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We show that the spectral function for single-particle excitations in a two-dimensional Fermi liquid 
has Lorentzian shape in the low energy limit. Landau quasi-particles have a uniquely defined 
spectral weight and a decay rate which is much smaller than the quasi-particle energy. By contrast, 
perturbation theory and the T-matrix approximation yield spurious deviations from Fermi liquid 
behavior, which are particularly pronounced for a linearized dispersion relation. 
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1. Introduction 

The low energy physics of interacting Fermi systems is 
usually governed by Landau's Fermi liquid theory, as long 
as no binding or symmetry-breaking occurs G|. Fermi 
liquid theory is based on the existence of quasi-particles, 
i.e. fermionic single-particle excitations with an energy- 
momentum relation £k that vanishes linearly as k ap- 
proaches the Fermi surface T of the system. In three- 
dimensional Fermi liquids quasi-particles decay with a 
rate 7k that vanishes quadratically near T ', as a conse- 
quence of the restricted phase space for low-energy scat- 
tering processes. The single-particle spectral function 
p(k, £) exhibits a Lorentzian-shaped peak of width 7k as 
a function of the energy variable £ for k close to T . 

In one-dimensional Fermi systems Landau's quasi- 
particles are unstable: the decay rate of single-particle 
excitations with sharp momenta is at least of the order 
of their energy and their spectral weight vanishes in the 
low-energy limit, giving rise to so-called Luttinger liquid 
behavior g|. 

No consensus has so far been reached on the proper- 
ties of single-particle excitations in two-dimensional Fermi 
systems. The anomalous properties of electrons moving 
in the CuCVplanes of high-T c superconductors have mo- 
tivated numerous speculations on a possible breakdown 
of Fermi liquid theory in 2D systems, at least for suffi- 
ciently strong interaction strength [|| and, according to 
Anderson Q, even for weak coupling. 

The development of new powerful renormalization 
group techniques J5], ^) has recently opened the way to- 
wards a rigorous non-perturbative control of interacting 
Fermi systems for sufficiently small yet finite coupling 
strength M. Significant rigorous results have already 



been derived for 2D systems. For example, the exis- 
tence of a finite jump in the momentum distribution has 
been established for two-dimensional Fermi systems with 
a non-symmetry-broken ground state JsJ . This excludes 
Luttinger liquid behavior in weakly interacting 2D sys- 
tems. 

So far, however, no rigorous results could be obtained 
for real-time (or frequency) dynamic quantities such as 
spectral functions. In particular, even at weak coupling 
the shape of the spectral function for single-particle ex- 
citations and the existence of well-defined quasi-particles 
are still under debate. Indeed, within second order per- 
turbation theory the imaginary part of the self-energy 
S(k, £) exhibits a sharp peak near £ = £k in two di- 
mensions ||, to be contrasted with the simple quadratic 
and k- independent energy-dependence of ImS(k,() in 
3D. Such a peak in the self-energy leads to a spectral 
function with two separate maxima instead of a single 
quasi-particle peak |9|, |l0| . However, analyzing an ana- 
lytic continuation of one-dimensional Luttinger liquids to 
higher dimensionality d, Castellani, Di Castro and one 
of us H have pointed out that such peaks in the per- 
turbative self-energy are smeared out in a random phase 
approximation (RPA) and probably also in an exact so- 
lution. Most recently a breakdown of Fermi liquid theory 
has been inferred from a divergence of the slope of the 
self-energy, computed within a T-matrix approximation 



(TMA), at£ = & y 



In this article we present explicit results for the self- 
energy computed within perturbation theory, TMA and 
RPA for a two-dimensional prototypical model with local 
interactions. We will then provide several simple argu- 
ments showing that the RPA, which yields a Fermi liq- 
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Figure 1: Second order Feynman diagram contributing to 

uid type result for self-energy and spectral function, is 
qualitatively correct, while unrenormalized perturbation 
theory and the T-matrix approximation produce artificial 
singularities in 2D systems. 



2. Perturbation theory and T-matrix 

We consider an isotropic continuum model with a local 
coupling as a prototype for Fermi systems with short- 
range repulsive interactions, 
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where a^ a and are the usual creation and annihilation 
operators for spin-^ fermions, £k is an (isotropic) disper- 
sion relation, j>0a coupling constant, and V the volume 
of the system. Local interactions between particles with 
parallel spin do not contribute due to the Pauli exclusion 
principle. A cutoff must be imposed for large momenta 
to make the model well-defined. 

The spectral function for single-particle excitations 
p(k, £) is directly related to the imaginary part of the 
one-particle Green function at real frequencies [jl2| 

G(k,0 
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Here /i is the chemical potential (fixing the particle den- 
sity). Note that /?(k, £) can be measured by angular re- 
solved photoemission. 

To first order in the coupling constant g, the self-energy 
is a constant which can be absorbed by shifting the chem- 
ical potential. Non-trivial dynamics enters at second or- 
der, where for a local interaction only a single Feynman 
diagram (see Fig. [I]) yields an energy-dependent contribu- 
tion £( 2 ) (k, £). The imaginary part, Im£( 2 ) (k, £), is cutoff 
independent and can be expressed as a two-dimensional 
integral, which can be easily performed numerically. Low- 
energy results for a quadratic (free-particle) dispersion 
£k = k 2 /2m and a linear dispersion £k = Vp|k| are shown 
in Fig. ^. The Fermi momentum kp and the Fermi ve- 
locity vp has been set one in both cases (vp = kp/m 
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Figure 2: Numerical results for Im£( 2 ) (k, £) at fixed |k| = 
l.Olkp for (a) a quadratic and (b) a linear dispersion 
relation, compared to (c) the asymptotic analytic result 
(||) . The fine-structure of the peak near £ = £ k is shown 
in the inset. 

for quadratic dispersion), and k has been fixed slightly 
outside the Fermi surface, i.e. |k| = l.Olkp. Note that 
the two curves almost coincide except very close to the 
point £ = £k = £k — A*, while the structure of the peak 
near £k depends sensitively on the curvature of £k near 
kp. Indeed, a strikingly simple analytic expression can be 
derived for the asymptotic low-energy behavior of ImS^ 2 ' 
(see Appendix): 
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for small £ at fixed k — £k/£ 7^ 1- F° r generic (non- 
linear) dispersion relations there are non-universal devia- 
tions from the leading asymptotic result in a tiny region 
of width around £ = £k- For k = hp (i.e. k = 0) one 
recovers the known perturbative result Im£( 2 ) (k^, £) cx 
£ 2 log|£| + 0(|£| 2 ) @. The logarithmic singularity for 
? ~ £k has already been noted earlier in the literature 
[^[ 10 1, but the complete asymptotic low-energy expres- 
sion (H) for the second order self-energy has not been 
found so far. 

As pointed out in another context by Castellani et al. 
H and also by Fukuyama and Ogata 0, a peak in the 
self-energy as in Fig. || would lead to a spectral function 
with two maxima centered around the quasi-particle en- 
ergy £k, which is reminiscent of the two-peak structure 
of spectral functions in one-dimensional Luttinger liquids 
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Figure 3: Numerical results for Im£(k, £) at fixed |k| = 
l.OOffci? within second order perturbation theory, TMA 
and RPA, for a coupling g = 1^ and a quadratic ek; the 
TMA requires a cutoff for large momenta which has been 
chosen as A = 5fcp. 

Several authors have computed the self-energy for two- 
dimensional Fermi systems within the so-called T-matrix 
approximation (TMA), i.e. summing all particle-particle 
ladder diagrams |Q , |f6| , p"7| . This approximation is ex- 
pected to be asymptotically exact at low density n — > 
in dimensions d > 2, where corrections are smaller by 
a factor n in 3D but only by a factor l/|log(n)| in 2D 
JTr| . Numerical results for the imaginary part of the T- 
matrix self-energy E T (k, £) , calculated for the model (|l]) , 
are shown in Fig. |3[ The peak at £ = £k is now bounded 
even for a linear ek- However, Yokoyama and Fukuyama 
[TT| have recently found that the real part of S r (k, £) has 
infinite slope at £ = £k for a linearized dispersion relation. 
Evaluating the renormalization factor 



z(k,o = [i-9£(k,0Mr 
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in the special limit k — > k^?, £ — > with £/£k ~> 1 they 
obtained the result " Z — > 0" , which was interpreted as a 
confirmation of Anderson's Q ideas about a breakdown 
of Fermi liquid theory in two dimensions. This interpreta- 
tion seems to be misleading. Firstly one cannot replace 
the function Z(k, £) by a single number " Z" if Z(k, £) 
has no (unique) limit for k — > k^, £ — > 0. Indeed, Z(k, £) 
vanishes o?iii/ in the above special limit. In particular, 
the jump in the momentum distribution function across 
the Fermi surface, which is also given by "Z" in a con- 
ventional Fermi liquid, does not vanish within the TMA 
in 2D systems! Secondly, one cannot trust a singular re- 
sult for a correlation function obtained from a truncated 
asymptotic expansion with respect to a "small" parame- 
ter, since higher order corrections are generally not uni- 
formly small for all momenta and energies. One should 



be even more worried in a case where corrections are sup- 
pressed at best by logarithmic factors 1/| log(n)| Jig ]. In- 
deed, the artificial sensitivity of the peaks in E(k, £) on ir- 
relevant non-linear terms in the dispersion relation found 
in perturbation theory and TMA indicates that the exact 
self-energy will look differently. 



3. Random phase approximation 

A self-energy with a regular energy-dependence which 
does not depend sensitively on irrelevant terms in ek is 
obtained from the random phase approximation 

E**(k,0=i j-j^Lj ^(q,.)G( >(k-q,e-c), 

(5) 

where is the non-interacting Green function and D 
the RPA effective interaction between particles with par- 
allel spin projection. For the locally interacting model, 
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where IT ' is the non-interacting polarization function. 
Numerical results for lvaE RPA are shown in Fig. 9. In 
contrast to the single sharp peak near £ = £k obtained 
in perturbation theory and TMA, the RPA result ex- 
hibits two smooth maxima with a width that is of the 
same order of magnitude as their height. These max- 
ima are due to long- wavelength collective charge and spin 
density fluctuations which contribute significantly to the 
RPA effective interaction; such collective fluctuations are 
not described by perturbation theory or TMA. The cor- 
responding spectral function has Lorentzian shape, with 
a smooth background that vanishes quickly in the low en- 
ergy limit. The quasi-particle decay rate 7k can thus be 
computed unambiguously from Im£ fl (k, £k) and turns 
out to be proportional to {k—kp) 2 log \ k— kp\ for arbitrary 
£k (while the perturbative ImS^ 2 ^(k, £k) is infinite for a 
linear ek, and finite else). The renormalization function 
Z(k, £) computed from the RPA self-energy has a unique 
limit Z for k -> kp, £ -> 0. 

The peaks in E(k, £) near £ = £k are produced by vir- 
tual excitations where the excited particles and holes have 
low energies and momenta near the Fermi surface close to 
the quasi-particle momentum k. In 3D systems the phase 
space for these particular excitations is very small and the 
low-energy behavior of ImE is dominated by low-energy 
particle-hole excitations with momenta all over the Fermi 
surface. Second order perturbation theory, TMA and 
RPA all yield the same simple behavior ImS(k, £) cx £ 2 
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with a possible weak k-dependence in the prefactor. In 
fact a renormalized perturbation theory where the bare 
interaction is replaced by the exact quasi-particle scatter- 
ing amplitudes yields the exact leading low-energy terms 
of ImS in 3D already at second order (in the scattering 
amplitudes) (l9) . For dimensionality d < 2 this is not true 
because virtual excitations with small momentum trans- 
fers dominate ImS Their contribution depends on 
the strength of the interaction vertex for small momen- 
tum transfers q (i.e. in the forward scattering channel) 
which is itself strongly renormalized by low-energy exci- 
tations Since the polarization function 11^ (q, uj) has 
contributions exclusively from low-energy states for small 
q, it is reasonable to reorganize the perturbation series 
by replacing bare interactions with RPA effective inter- 
actions. The simplest approximation for the self-energy 
is then to dress the bare fermionic propagator by only 
one effective interaction. It has been pointed out already 
earlier by Eschrig et al. (2(]] that all RPA diagrams yield 
equally important contributions to S in a renormalized 
low-energy expansion in two dimensions. 

The random phase approximation is not an arbitrary 
resummation of Feynman diagrams, but is controlled by 
the expansion parameter 1/N, where N is the number 
of spin states in a generalization of model (|l|) , where in- 
stead of two possible spin-projections one allows for N 
different spin states. Rescaling the coupling constant by 
a factor 1/N, one can take the large N limit and classify 
Feynman diagrams by powers of 1/N. Since each interac- 
tion line in a diagram yields a factor 1 /N and each closed 
loop a factor N, all non-RPA contributions are of order 
(1/N) 2 . In contrast to plain perturbation theory and the 
TMA, which are controlled by the coupling constant g 
and factors 1/| log(n)|, respectively, the RPA yields a reg- 
ular result for £(k,£). Thus 1/N seems to be a suitable 
expansion parameter for 2D Fermi liquids, with higher 
order corrections being uniformly small for all k, £ (for 
large N). 

Even more convincing is a comparison with one- 
dimensional Fermi systems, where exact results are avail- 
able 0] . Since the peaks in the self-energy discussed above 
are due to excitations of particles and holes with almost 
equal momenta, it is sufficient to consider a special case 
of the exactly soluble Luttinger model where only the 
particles and holes near the right (or left) Fermi point 
interact. In Fig. |] we show the exact result for ImS(fc, £) 
for this model, compared to the RPA result and the result 
of second order perturbation theory. The sharp peak (a 
(5-function for a linearized dispersion) appearing in per- 
turbation theory is replaced by a smooth maximum with 
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Figure 4: ImE(fc,£) at fixed k = 1.001/cf for a one- 
dimensional model with a linear dispersion relation and 
a coupling g — Ivf] the exact result is compared to per- 
turbation theory and RPA (the vertical arrow represents 
a 5-function). 

equal height and width in the exact solution. The RPA 
correctly captures the width and height of the exact func- 
tion, though not its form. One may thus guess that the 
two broad maxima in ImE(k, £) obtained from the RPA 
in 2D have to be replaced by one broad maximum in the 
full solution. In any case the spectral function will have 
Lorentzian shape. Note that the RPA correctly signals 
all the non-Fermi liquid singularities in one dimension 
(lmS RPA (k, £) is linear in k—kF, £ in ID), while it yields 
regular Fermi liquid behavior in 2D. 

Our explicit (numerical and analytical) results have 
been obtained for an isotropic model with a weak local 
interaction. The qualitative conclusions hold more gen- 
erally for any Fermi liquid with short-range interactions 
and a Fermi surface with finite curvature, as long as no 
instability associated with diverging renomalized interac- 
tions occurs. 



4. Conclusion 

In summary, we have shown that the spectral function for 
single-particle excitations in a two-dimensional Fermi liq- 
uid with short-range interactions has Lorentzian shape. 
Quasi-particles have a uniquely defined spectral weight 
Z (equal to the jump of the momentum distribution 
across the Fermi surface) and a decay rate proportional 
to k' 2 logfc', where k' is the distance of the quasi-particle 
momentum from the Fermi surface. The latter result is 
not new, but is has so far been derived only by "good 
luck", i.e. by computing Irn£(k, £k) within perturbation 
theory (2l| or TMA |15| for a quadratic dispersion, where 
the spurious peak in ImE(k, £) reaches its maximum at an 
energy of order ^/cf below £k! We have provided several 
arguments showing that the random phase approximation 
yields a qualitatively correct result for the spectral func- 
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tion, while plain perturbation theory and the T-matrix 
approximation produce spurious singularities. 
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Appendix: 

Self-energy in second order perturbation theory 

For the continuum model ([!]) the second order 
contribution to the (time-ordered) self-energy can 
be calculated using the two-particle propagator 
K(p,u) = iJ$j s fg G(°) (p - k, u - G(°) (k, . 
In the following we make use of the spectral represen- 
tations of the self-energy, E(k,£) = J d£'S(k, £')/(£ - 
f + iO + sgn£'), the free one-particle propagator, 
(k, f ) = / ' p (0) (k, £')/(£-£' + iO+sgn £'), and the 
two-particle propagator, K(p, cj) = J" dw' A pp (p, a/)/ (cj — 
u/ + iO + sgno/). For a quadratic dispersion £k = k 2 /2m, 
App is given by pi 
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where cjo := p 2 /4— 1, cj± := p 2 /2±p, and the free density 
of states Np — kp/2Trvp. Here and in the following mo- 
mentum variables are measured in units of kp and energy 
variables in units of vpkp. 

The second order contribution to the spectral function 
of the self-energy can be written as 

ftduj f d 2 p 
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is the angle-integrated imaginary part of the free propa- 
gator G^ ) , <f> being the angle between p and k |l(| . 

For small £ <C 1 the integration boundaries \k — 
V2C + 1| < p < k+^/2( + 1, given by the step- function 6 
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Figure 5: Integration boundaries in the (p,u>) plane for 
< £ <C 1 and e := k — kp > 0. p, e in units of kp ; w, £ 
in units of vpkp. 



in eq. (A.3), can be linearized to \k — £ — 1| < p < fe+C+1. 
The same can be done with the boundaries coo, u>± in eq. 
(A.l). This leaves us with trapezoid integration regions 
as given by the shaded areas in Fig. Hfor the case £ > 
and k — kp > 0. Also shown in Fig. |5| are subareas (1) 
- (3), which can be defined analogously for other choices 
of £ and k. In these subareas we can use the following 
asymptotic expressions for p and A pp (e := k — kp): 
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(A.4) 

if (p,C)e(io) 
if (p,C)G (1), (16) 
if (P,C)e(2) 

if (p,C)G (3). 

(A.5) 

These approximations lead to a result for the self-energy 
which is exact within second order in £ at fixed n := 

we * 1- 

With these asymptotic expressions the integration over 
p can be performed analytically. The results for the dif- 
ferent subareas are rather lengthy, but can be simplified 
by a further expansion in u> = 0(£) and e = 0(£), keeping 
the ratio k = e/u> fixed. Now also the integration over u> 
can be performed analytically, resulting in the simple ex- 
pression eq. (|) (Note that S(k,£) = |Im£(k, £)|/tt). The 
main contributions to S come from the momenta p ~ 
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and p w 2k p . The contributions from p sa are 
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causing the logarithmic divergence at k = 1. 
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